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ABSTRACT: 

In the framework of heavy baryon chiral perturbation theory, we give the chiral ex¬ 
pansion for the ttA —> tittN threshold amplitudes Di and D 2 to quadratic order in the 
pion mass. The theoretical results agree within one standard deviation with the em¬ 
pirical values. We also derive a relation between the two threshold amplitudes of the 
reaction nN —>■ titiN and the tttt S-wave scattering lengths, Og and Og, respectively, to 
order 0{M^). We show that there are uncertainties mostly related to resonance exci¬ 
tation which make an accurate determination of the tttt scattering length Og from the 
TTTrA threshold amplitudes at present very difficult. The situation is different in the tttt 
isospin two final state. Here, the chiral series converges and one finds Og = —0.031 ±0.007 
consistent with the one-loop chiral perturbation theory prediction. 





1 Introduction 


Elastic pion-pion scattering in the threshold region is the purest process to test our 
understanding of the chiral symmetry breaking in QCD. Already in the early days of 
current algebra, Weinberg [|l| showed that the tttt S-wave scattering lengths Qq (with 
1 = 0,2 the total isospin of the two-pion system) vanish in the chiral limit of zero quark 
masses. In particular, he predicted Og = 7Ml/{?>2 'kF^) = 0.16 and Og = (—2/7)ao. This 
prediction was further sharpened by Gasser and Leutwyler [Q in the framework of 
chiral perturbation theory (CHPT), which is the effective held theory of the standard 
model at low energies. In Ref.|p a very accurate prediction for the isospin zero, S-wave 
scattering length was given, Og = 0.20 ±0.01, which amounts to a 25% increase compared 
to the current algebra value. The main assumption underlying this result is that the order 
parameter of the chiral symmetry breaking, B = — < 0|gg|0 > /Fl (with < 0|gg|0 > 
the scalar quark condensate and ~ 93 MeV the pion decay constant) is considerably 
larger than i%, B 3> which follows e.g. from the standard analysis to determine the 
light quark mass ratios from the Goldstone boson masses. Another scenario, in which the 
quark condensate is very much smaller and consequently i? ~ has been discussed in 
Refs. 0 [§] . In this approach, the ratio of the average u and d mass to the strange quark 
mass is decreased, typically 2ms /± fnd) < 10 and the resulting scattering length Og 
increases, typically Og > 0.27. To settle this very important issue, it is mandatory to 
determine the tttt S-wave scattering lengths within an accuracy of about 20% (or better). 
For a review on these topics, see e.g. p. 

It is, however, not straightforward to determine the Trvr phase shifts in the threshold 
region experimentally. A few possible candidates are Ki^-decajs, pionic molecules or the 
reaction ttN ttitN. It is this latter process we will be concerned with in the following. 

To be more precise, consider single pion production in the threshold region and above. 
Already Weinberg pointed out that the one-pion exchange diagram contains the four- 
pion vertex (with one pion leg off-shell). This opens two possibilities of extracting the 
on-shell tttt interaction. First, one can consider peripheral processes at higher energies 
but low momentum transfer and try to isolate the pion-pole by standard Ghew-Low type 
techniques (for recent work in this direction, see e.g. |0). Here, we will concentrate on the 
second way, namely to directly relate the two independent titiN threshold amplitudes to 
the scattering lengths Og and Og. This approach was pioneered by Olsson and Turner 
and has been used ever since in most analyses of the threshold iiN ttitN data. However, 
the Olsson-Turner approach predates QGD, it is not applicable any more beyond tree 
level. In particular, in its original formulation a parameter which is a measure of the 
type of chiral symmetry breaking, is left free. In QGD, this parameter ^ is exactly zero. 
A critical discussion of these topics can be found in . 

On the other hand, over the last few years an impressive series of experiments have 
measured the total cross section for the processes ttN —>■ tittN quite close to threshold |TT| , 
0 1^, [^, 1^. Extracted values for the tttt scattering lengths are based on the Olsson- 
Turner approach with ^ ^ d [O, M. Therefore, it is necessary to work out a more 


precise relation between the threshold tittN amplitudes and the tttt S-wave scattering 
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lengths beyond tree level. A first step in this direction was made in Refs.[17| |18] where 


an improved low-energy representation accounting for corrections of order to the tree 
level relations was formulated. This led to novel low-energy theorems, which can be 
directly compared with the threshold data. Not unexpectedly, one hnds a satisfactory 
description for the channel 7r~^p 7T~^7r~^n and some signihcant deviations for the process 

7r~p —>■ TT^TT^n, reflecting the relative weakness/strength of the pion-pion interaction in a 
state with isospin two/zero (in the S-wave). Here, we present the results of a calculation 
to one loop accuracy in heavy baryon chiral perturbation theory (HBCHPT) p0| , 
which accounts for all corrections up-to-and-including 0{M^) to the tree level result 
and which is therefore sensitive to the one-loop corrections to the tttt scattering lengths, 
besides many other contributions. 

The paper is organized as follows. In sections 1,1 and § we formulate the problem, 
give necessary kinematics and discuss briefly the effective Lagrangian that will be used. 
Section ^ contains the principal results of this paper, namely the Born, one loop and 
counterterm contributions to the two threshold tittN amplitudes up-to-and including 
order M^. Numerical results are discussed in section ^ and a short summary is given in 
section Some technicalities are relegated to the appendices. 


2 Prelude I: Threshold kinematics for ttAT —> ttttAT 


Consider the process 

7r“(A;) + N{pi) 7i\qi) + 7i%q2) + N{p2) , (1) 

with 'a, b, d pion isospin indices. N denotes the nucleon (neutron or proton). At threshold 
and in the centre-of-mass frame, we have qi = q 2 = {Mn, 0, 0, 0), with the pion mass. 
Using the pseudoscalar quark density P“(x) = q[x)i'')^T^q{x) as interpolating pion held, 
standard LSZ reduction leads to 

- J d^xd^y < N\T[P^{0)P\x)P%y)]\N > 


with dehned via 




(M2 - ql) (M2 - ql) (M^ - k^) 


T 


off-shell 


< 0|P“(0)|7r'’ >= ^“''G 


( 2 ) 

(3) 


At threshold, the on-shell amplitude in the 'k°‘N centre-of-mass system can be expressed 
in terms of two threshold amplitudes. 


t: 


on—shell _ 


cms 


= T = id ■ k 


Pi(r'’(5“"+ Pr'’) + 


rfec 


(4) 


The quantities D 12 in eq.(^ are related to the commonly used amplitudes A 2 ij^^, with 
I the total isospin of the initial nN system and the isospin of the two-pion system in 
the hnal state, via 


A32 — \/T 0 Pi, ^10 — — 2 Di — 3D2 


(5) 
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What we are after is the chiral expansion of the Di and D 2 . It takes the form (with V a 
generic symbol for D 12 ) 

V = fo + fifi + + ■■■, 11 = M-^/m , (6) 

modnlo logarithms and we have introdnced the pion to nncleon mass (m) ratio, /r. 


3 Prelude II: Evolution and formulation of the prob¬ 
lem 


The modnlns of the threshold amplitnde is determined by the extrapolation of the mea- 
snred total cross section in the threshold region via 


\A{TnTN)\^ 


a{nN —vrvrV) 


(7) 


where is the incident laboratory pion kinetic energy, S' is a Bose symmetry factor 
(S' = 1/2 if the final two pions are identical, otherwise it is nnity), and 


C=Ml (^j^)y3(2 + ^)>/y2 + 3,j)‘/yi + 2f,)-“''V (8) 


The threshold modnlns has been obtained in this way for the five charge states initiated by 
TT^p 1^. Explicit isospin violation dne to the electromagnetic mass differences has been 
removed throngh the kinematics of the threshold valne and the threshold amplitnde 
modnlns is assnmed to be isospin invariantBy Watson’s theorem the threshold 
amplitude has the phase of the initial elastic P-wave vriV scattering amplitude 

(up to an overall sign). The threshold production amplitude complex phase is then ~ 
—4° for initial vriV isospin 3/2 and ~ 2° for initial isospin 1/2, as given by the respective 
phases evaluated at the cms momentum of 213.6 MeV (the ttttV threshold). The threshold 
production amplitudes are thus nearly real. At threshold the final tttt S-wave state must 
have isospin 0 or 2 by extended Bose symmetry and hence there are only two independent 
threshold amplitudes as discussed before. From the measured process amplitude moduli 
a unique value of ^10 and Az 2 can be found up to an overall sign. We note that the sign 
is fixed by the chiral expansion (as discussed below). 

In the Olsson-Turner approach [Q (with .^ = 0), it follows that 


^32 

klio 


Q-kN ' ^0 

-MI 


-2^/W7^— ^ + d2 


dvr 


QwN 

m 


m 


(9) 


with gfjrTv = 13.4 the strong pion-nucleon coupling constant. The above result is a con¬ 
sequence of the dominance of the pion exchange and contact diagrams. To lowest order, 

more systematic study of isospin violation is certainly needed. 
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the two “shift” constants h/ arise from the ” sub-leading” diagrams involving three pion 
absorptions/emissions on the nucleon line, compare Fig. 1. The dj are of order In 

the context of QCD, the relations Eq.@ are equivalent to the tree level CHPT results if 
the ai are the tree level tttt scattering lengths a la Weinberg. The improved representation 


of [0, which is based on the hrst corrections to the Olsson-Turner result, takes the 


form 


*^32 — 



m 

AlO = 

(1 + 
m 


(1+2/^) 
37 
14 ^ 






( 10 ) 


where the new shift constants dog have the form 


d/ — dq -|- d\MT^ d'jM^ , 


/ = 0,2 


( 11 ) 


modulo logs. One notices that the correction of order M^- is comparable in size to the 
leading term (approximately 40% and 50% for ^10 and ^ 1 . 32 , respectively). Therefore, it 
is mandatory to calculate (at least) the coefficients d?. Also, at that order the one-loop 
corrections to the S-wave tttt scattering lengths appear ^ The problem investigated 
here is thus nothing but the calculation of the constants d? in eq. 


4 Prelude III: Effective Lagrangian 


In this section, we will briefly discuss the chiral effective pion-nucleon Lagrangian as well 
as the pionic one underlying our calculation. Many additional details are spelled out in 
Refs.0 To explore in a systematic fashion the consequences of spontaneous and 

explicit chiral symmetry breaking of QCD, we make use of baryon chiral perturbation 
theory (in the heavy mass formulation) (HBCHPT). The nucleons are considered as 
extremely heavy. This allows to decompose the nucleon Dirac spinor into ’’large” (id) 
and ’’small” {h) components 


^(x) = + h{x)} , 


( 12 ) 


with the nucleon four-velocity, = 1 , and the velocity eigenhelds are dehned via 

1 / H = H and 1 / h = — h.|^ Eliminating the ’’small” component held h (which generates 
1/m corrections), the leading order chiral vriV Lagrangian reads 


= H{iv -D+IaS- u)H . 

Here the pions are collected in a SU(2) matrix-valued held U{x) 


U(x) = - 


\JF'^ — Ti {xy + ir ■ 71 (x) 


(13) 


(14) 


#^The role of is to single out a particular reference frame 1^, here the 7r“A^ centre-of-mass frame. 
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with F the pion decay constant in the chiral limit and the so-called cr-model gauge has 
been chosen which is of particular convenience for our calculations in the nucleon sector. 
In eq. ([T^) = d^ + V^ denotes the nucleon chiral covariant derivative, is a covariant 

O 

generalization of the Pauli spin vector, 9a the nucleon axial vector coupling constant in 
the chiral limit, = iu^VfJJu^ with u = \/U and the covariant derivative acting on 
the pion helds. To leading order, 0{q), one has to calculate tree diagrams from 

+ ^Tr{V'^f/V;,17^ -I- x+} , X+ = + ux^u (15) 

where the second term is the lowest order mesonic chiral effective Lagrangian, the non¬ 
linear cr-model coupled to external sources. The quantity y contains the light quark 
mass m and external scalar and pseudoscalar sources (the latter are actually needed to 
compute correlators of the pseudoscalar quark density as in eq. (H) ) . Later, we will need 
the hve-pion-nucleon vertex. Expanding u in powers oi (j) = r ■ n/F gives 


u = 1 + -6 — 


8 


+ 


16 


128 


+ 


7i 


+ 0{^«) , 


(16) 


and Ufj, follows correspondingly, = i{u\dfj,u}. 

To one-loop accuracy, i.e. order 0{q^), one has to consider tree graphs from 


— 




I ^(2) I ^(3) I ^(2) I ^(4) 

'’^'TT'TT ^^'TT'TT 


’( 3 ) 


-ttW 


-ttN 


(17) 


i‘2') 

where the structure of is discussed in detail in 
,( 2 ) 


1 l 8 | and, on a pedagogical level, in 


All terms in are hnite. The hrst divergences appear to 0{q^) in HBCHPT, the 


corresponding determinant has been worked out by Ecker 


^(3,div) 


- 3- 

Y^hHOiH , 6, = &[(A) + ^L 

i=l ^ 


(18) 


where the Oi are monomials in the helds, A is the scale of dimensional regularization and 
the bi differ by a factor ( 47 rTV )2 from the ones in [^, and 


L = 


A 


d-i r 


IGtt^ 


- 7e + 1) 


(19) 


(3) 

with jE = 0.5772 the Euler-Mascheroni constant. There are also terms in with hnite 
coefficients. The corresponding low-energy constants will be estimated by resonance 
exchange. It is important to note that some of the terms in simply 1/m and 

l/m2 corrections from the original Dirac Lagrangian, like e.g. HD^/{2m)H (for details, 
see 1^). contains terms proportional to the low-energy constants Ci, 02 , 03 , 04 . The 
latter are related to the ttN cr-term and ttN scattering lengths as discussed below. In order 
to restore unitarity in a perturbative fashion, one has to include (pion) loop diagrams. 
In HBCHPT, there exists a strict one-to-one correspondence between the expansion of 
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any observable in small external momenta and quark masses and the expansion in the 
number of (pion) loops. In what follows we will work within the one-loop approximation 
corresponding to chiral power 0{q^). To obtain all contributions at order one has to 
supplement the chiral effective Lagrangian by the additional term [^. It serves to 
cancel some of the divergences of certain loop diagrams and contains the mesonic low- 
energy constants f'i, 2 , 3 , 4 . The latter encode information about the chiral corrections to the 
TTTT scattering lengths. We use the following form of |^ , 


= ^(TrV^f/V'^f/^)' + ^Tr(V^[/V.17^)Tr(V^f/V"17^) + Y|(TrX+)' 

-h^{2Tr(V^f/V'^f/t)Trx+ + 2Tr(xtf/;^t^ + ^^t^f/t)_4Tr(;^t^)_(Trx_)2| + ... (20) 

where the ellipsis stands for other terms of order q'^ which do, however, not contribute in 
our case. The hnite pieces of the low-energy constants in eq.(^0D are renormalization 
scale dependent and are related to the li of ref. [0] via 


4 = 967r2£^(A) -21n- 
4 = -647r2£^(A) -21n. 


4 = 487r2£^(A) -21n 
h = 167r2£^(A) -21n 




A 


A 

M, 

T ’ ■"* A 

and their actual values will be discussed later. From the Lagrangian eq. (pOD , one derives 
the chiral corrections to the S-wave tttt scattering lengths [§] (we do not exhibit the 
explicit scale-dependence of the any more) 


( 21 ) 


flO ^ ^ 

° 327rF2 

o Ml 


7 Ml r f M, 
1 + ' 


2/5 


9 In 




2M2 

+ ^ 


Oq — 


IQiiFl 


1 + 


( 




, /fain—--I + 

VdTTF^y V A 2/ 


7F2 
2Ml 


F2 


-Ml - All - ^3 + ^4 


( 22 ) 


(23) 


with Mjr and F^^ the empirical values (i.e. the corresponding chiral corrections have been 
accounted for, see also |^). We have now assembled all tools for calculating the one loop 
corrections for tiN —> imN. 


5 ttAT —> ttttAT to one loop 

Before presenting the results of the calculation, some general remarks are in order. The 
various contributions to the chiral expansion of the invariant functions Di ^2 to order q^ 
can be grouped as 

•p(3) _ pBorn _|_ ^^one-loop _|_ ^ct ^24^ 

where "pBom subsumes the lowest order relativistic tree graphs and all kinematical correc¬ 
tions to it (which are suppressed by powers of 1/m), the generic one loop graphs 
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and the counter terms, which absorb the divergences from the loops and there are, of 
course, additional hnite ones. At threshold, the calculation simplihes since we have the 
selection rules 

S ■ qi = S ■ q 2 = V ■ {qi - q 2 ) = ^ , v ■ qi = v ■ q 2 = v ■ k = 2M^ + 0{—) . (25) 

m 

Also, from the start we will perform mass and coupling constant renormalization, as 
explained in the context of the Born graphs (see below). We will give a fairly detailed 
description of the renormalization of the remaining divergences since that will serve as an 
excellent check on the calculations. Furthermore, the extraction of the various low-energy 
constants and a thorough discussion of the related uncertainties is mandatory to really 
hlter out the sensitivity of the threshold tittN amplitudes to the S-wave tttt scattering 
lengths. 


5.1 Renormalized Born terms 

From the three-pion-nucleon seagull and the pion-pole diagram]^ one immediately hnds 
the leading 0{q) contribution to Di ^2 


Di = 


9a 

8F3 


Do = 


^9a 

'8F3 


(26) 


It is then most economic to calculate the relativistic tree graphs and expand the result 


in powers of /i, see Fig. 1. This gives automatically all kinematical corrections to eq. (^6[) 
and reads to order 0{q^) 


DBovn ^ 


9a 

8F3 


7 

1 + 


1 3 2 

8 + 2®-* 




9a 

8F3 


-3- 


17 2f 

2-® + " ( 


11 


+m 


(27) 


In what follows, we have to renormalize the factor 9a / (i.e. the chiral limit value) to 
the physical value g^N/ifuF^). This renormalization procedure subsumes a host of loop 
and counter term corrections. For doing that, we hrst have to consider the pion mass, 
decay constant and so on 0. To one loop we have for the pseudoscalar coupling (cf. 
Fig. 2) 


G 


TT 


G 




1 

IbTT^ 



(28) 


with G the chiral limit value of Similarly, we hud for the pion decay constant 


F^ 



1 

8n‘^ 



(29) 


"^3 It is important to remark that this splitting has no physical meaning. The contributions of the 
threshold amplitudes D12 coming from the seagull and the pion pole depend on the choice of interpolating 
pion field, i.e. how one parametrizes U{x) through some pion field. The sum as a physical quantity is of 
course unique and independent of the particular choice. 
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and the pion mass renormalization reads 


Ml = 


mV 1 

1 + + 




(30) 


with = 2mB the leading term in the quark mass expansion of the pion mass squared. 
With that, the pion propagator takes the form (in the a-model gauge), 


i 


M2 


q2 _ My ^ F2 


2L + 24 + 


1 , M 
—^ In “w 
Svr^ A 


(31) 


The pertinent diagrams contributing to the renormalization of the pion-nucleon vertex 
are shown in Fig. 3. The appropriate one-loop graphs for vriV —> ttuN which account for 
this renormalization will have to be identihed and subtracted accordingly, see the next 
section. The coupling constant renormalization can be written as 


9nN _ ^4 _ ‘^biiMl 

m V- V FI 


(32) 


where the constant hn is hnite. Its value is hxed from the known Goldberger-Treiman 
discrepancy. However, there remains a hnite contribution to Di ^2 from the corresponding 
chiral power three Lagrangian, 


cfl = bnp, H iS ■ (Dx-) H 


after the Qt^n renormalization, eq.(^), has been performed, 

= 0, 

This has to be accounted for. 


nGTR _ ^ i. 

F>2 - -9A Oil • 


(33) 


(34) 


5.2 One loop graphs 

In Fig. 4, we show the 36 different one-loop diagrams that contribute at threshold (we do 
not display graphs in which the two out-going pions are interchanged, 6 c). We have 


made use of the selection rules, eq. (p5|) , and omitted all those diagrams which according 
to these rule are equal to zero. This means that if one wants to extend this calculation 
to kinematics above threshold, one has to consider many more diagrams than shown 
in Fig. 4 since then the selection rules do not apply anymore. Of course, many of the 
graphs shown contribute to mass and coupling constant renormalization. Concerning 
the chiral corrections to the pion-pion interaction, the interesting diagrams are the ones 
numbered 16, 17 and 18. The calculation of all these diagrams is straightforward but 
somewhat tedious. It is most economically done in the basis of the loop functions dehned 
in appendix B of ref. [Q. In appendix A, we assemble some novel loop functions not 


considered in 18 . 
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Putting all pieces together, one has (after renormalization of and 


£^loop _ 9A-^^ 


2 r 


327r2F5 


9\ 


In 


A 


24 6 ^2 


9 aMI \g\ 1 


^loop 


9aMI r,13_i(, 2),,, 

4 A 24 


327r2F5 


F5 L2 
73 , 16(71 


^)v/3 ln (2 + Vs) - 4 + 10 ^ 

o 4 


(35) 


1 ^TT 

■(%A-2)^^’^(2 + A/3)+iY^5'A + 4/ 


FI 




with 


r-l 


X 


1=1 djX- 

'0 G(1-x)(1 + 3x) 


arctan 


X 


(1 — x)(l + 3x) 


= 0.6456 


(36) 

(37) 


The imaginary part in eqs.(^,^) is due to the diagrams numbered 20, 21, 24,..., 28. For 
these pion-nucleon rescattering type of graphs, the pertinent loop functions have to be 
evaluated at a; = 2 M 7 r, which is well above the branch point cuo = M^r. This is similar 
to the effect observed in the calculation of the threshold amplitudes for the reaction 
yiV —>■ tittN, where one also hnds an imaginary part at threshold |^. If one calculates 
from eqs. (^,^) the imaginary parts of the isospin amplitudes Alio and AI 32 and compares 
them to those demanded by Watson’s theorem (using the experimental £t values for the 
real parts), one makes the following observation. The phase is approximately correct 
for the isospin 3/2 case but an order of magnitude too large with the wrong sign in the 
isospin 1/2 case. The reason for this is that the tree level phases which we encounter 
here satisfy = 4^31 which does not hold for the empirical P 2 /,i phase at the vrTriV 
threshold. Nevertheless, the appearance of the tree level nN phases in eqs. ( |35| , |36[) serves 
as a good check on the one-loop calculation. We encounter here the standard problem in 
chiral perturbation theory that for getting the imaginary parts better, one has to perform 
a higher order calculation. As discussed in the introduction, D 12 have to be almost real, 
therefore we will neglect in the following their imaginary parts. Next, we have to perform 
the remaining renormalizations to get rid of the terms proportional to L. 


5.3 Renormalization 

We proceed in two steps. First, we consider the divergences related to the pion-nucleon 


('Q\ 

Lagrangian eq.([18D. The following operators as dehned in Ref.[^] give a non 


vanishing contribution at threshold; O4, O5, Oq, O 7 , Og, On, Ois and O 2 o- Some of these 
are equivalent at threshold, these are O 4 and Oq, O 5 and O 7 and the combination of 
Oi 7 plus Ois- The corresponding /3’s are = 9 Ai)- — 9 ]dl‘^-, Pi + = — 9 aI‘^ 

and Pn + P 18 = (2 — 3(7^)/4. The operator Og = S ■ uTt{x+) h^is two terms with three 
pion helds. First, there is a three-pion vertex from u^. This contribution is, however, 
completely contained in the renormalization of Second, there is a term with one 
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pion from and two from Tr(x+) with /Sg = (7^(4 — g\)/ 8 . Finally, there is the operator 
O 20 = iv ■ DTr(x_|_) + h.c.. The relevant contribution has a vertex with two pious coming 
from Tr(x_|_) with j 32 Q = —9g\/lQ. Adding up all these counterterm countributions which 
in the case of O 17 , Ois and O 20 come from two step processes with a. a ■ k interaction for 
the incoming pion 7r“, we have (the scale-dependence of the is not made explicit) 


£)(ct,3) 




&5 + &7 + 


gA 

2 F2 

TT 


(1 - gl) L 


(38) 


Diet,3) 


Ml 


2(6; + bl + bl) - 49.4(6;, + b\, + 6;„) + ||( 5 gJ - 2) L 


(39) 


where the terms ^ g\L are cancelled by the inhnities in and - 02 °°^, 

respectively. So we are left with the following divergent pieces: 


eqs. (^p 6 D , 




div,loop+ct3 _gAMl.?, 


F5 


L) 


^div,loop+ct3 _ 


gAM^ 27 

P5 ^ 8 


L) 


(40) 


which have to be cancelled by the counter terms from ClJ. After renormaliztion of F,,, 
Gt, and gT,N, the total contribution from reads 


^{ct,4) 




D 


(ct,4) 

2 


gAH 

F5 


2 r 





(41) 

(42) 


We remark that the low-energy constant £3 does only appear via the renormalization of 
the pion mass with the appropriate insertion from CH for the pion hooking on to the 
nucleon (cf. graph 12 in Fig.4) since [Tr(x+)]^ has no four-pion vertex in the a-model 
gaugeComparison of eqs. (|4l|j4^) with eq.(|40|) leads to the desired result, namely 


Df" = + Df"-* = 0 , ! = 1,2 . 


(43) 


This cancellation of divergences serves as an important check on our calculation, 
hnite counter term contribution can be compactly written as 


ct,fin _ M^ 


^^ct,nn _ 


F3 

-TT •- 


iA 

F 2 


()r pr o 

A — A. — Zf 




^2 


K 

F3 


2 r 


U 

>2 


+ + +s;. 


The 


(44) 

(45) 


where the S[ 2 subsume the appropriate contributions of the plus additional hnite pieces 
from cH and from 1/m suppressed corrections from cH- The estimation of these hnite 
pieces will be discussed in the next section. 

have checked that in other parametrizations of U where [Tr(y+)]2 has a four-pion vertex, the 
final result is the same. 
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5.4 Finite contributions and estimation of low—energy constants 


The most difficult task is to pin down the finite terms 61 2 iii egs-f^^Sf). These can be 
split into two distinct contributions. First, there are 1/m suppressed terms with insertions 
from the relativistic chiral order two Lagrangian In Ref.||17|| we had already shown 
that such terms cancel at order which allowed to formulate low-energy theorems for 
Di 2 independent of the corresponding low-energy constants q. Here, we are working 
one order further and thus such contributions appear, some examples are shown in Fig.5. 
These are operators of dimension three which contribute at threshold. The other type 
of terms are related to the values of the various &[(A) due to the renormalization and 
additional finite ones from In the absence of a complete data set to fit these, we 

will make use of the resonance saturation principle. This procedure will induce some 
uncertainty in our final results, see the discussions in [Q and 
this is the term ~ bu discussed in section ^ 
explained by resonance exchange. 

Consider first the contributions due to insertions from These can be either calcu- 


The exception to 
This particular contribution can not be 


lated by performing the path integral as in ref. and expanding to the desired order, or, 
more economically, by using the relativistic Lagrangian It is this latter method we 

IPS] 


are using.|3 It is same method used in the calculation of the 1/m suppressed (kinemat- 


ical) corrections to the tree graphs, compare section |5.1| . It is based on the observation 
that all operators which have fixed coefficients like 1/m, 1/m^, and so on are noth¬ 

ing but the expansion coefficients of the relativistic theory (see e.g. Refs. p0[| , p3| ). The 
corresponding dimension two effective relativistic pion-nucleon Lagrangian reads 10 0 


c 


(2,rel) 

tzN 


= ciTTTr(x+) Du 

4m 




^ D^,Du ^Tr(n^M^) 




(46) 


where T denotes the relativistic nucleon field and the ellipsis stands for other terms not 
needed here. The c* are normalized such that we can identifiy them with the correspond¬ 
ing low-energy constants of the heavy baryon Lagrangian (truncated at order q^) (for 
definitions, see e.g. |]T^). Note that the constants C 2 and C 2 are related to the C 4 and C 5 
of App. A in [^] and that they have been renamed in comparison to Ref. ||27||. To leading 


order in the heavy baryon Lagrangian, we have C 2 = C 2 -|- C 2 in the notation of Ref. 
The contribution of these terms to the Di ^2 follows as. 


^(cp _ 9aMI 


mF^ 


D 
2 


-2,Ci T 7^C2 -|- 2c2 “1“ 2c3 T ~C4 


1 

—( 
2 


(47) 


= 


9aM^ 

mF^ 


C2 ~h 2 c2 C4 


(48) 


appendix we derive the result for D12 directly from the path-integral formulation of the 
heavy-nucleon Lagrangian. 
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For the numerical evaluation of the Cj and to the accuracy we are working, we can fix 
them to order (i.e. just calculating tree graphs) from available pion-nucleon scattering 
data. Therefore, the values presented here will differ from previous estimates which 
included contributions. This difference is, of course, just one of the many corrections of 
0 {M^) in Di and D 2 and only becomes important if one wants to extend the calculation 
presented here to the next order. First, we invoke the subthreshold expansion of the 
standard invariant pion-nucleon amplitudes with the pseudovector Born terms subtracted 
(as indicated by the ’bar’) 

X = E Xm,n r , X = {kl+, 5+, A-,B-} , (49) 

m^n 


with t the invariant momentum transfer squared, v = {s — u)/4:m {s,t,u are the con¬ 
ventional Mandelstam variables subject to the constraint s + t + u = 2M^ + 2m?) and 
fc = 1 (0) if the function considered is odd (even) in z/. Retaining terms to order (9(z/^,t), 
one finds from eq.( 


R+ = 


4ciM^ C3 ,0 \ 2 c2 2 


2r' 

' p 2 


A- = 


2 mc4 

F? 


,B- = 


+ 


2 c4m 


2p2 p2 


(60) 


where the first term in B stems from the celebrated Weinberg-Tomozawa term 
From that, we deduce the following relations to order g^. 


31 


®oo ~ 


2 Ml 

TT 

~pF 


(c3-2ci) ,aoo = 


2 mc 4 

F 2 


5 ®10 ~ 


F 2 


5 ®01 ~ 


C3_ 

F 2 


~ p2 ’ ^00 ~ 2^2 4c4m) . 


Further information is obtained for Ci from the irN-a term [^, a^Ar(O) = 
(to order g^) and for C 3 and C 4 from the wave scattering volumina Cq and di, 


Co = 


C 3 


27iF? ’ 


di = 


47rF2 


C 4 4- 




AmJ 


(51) 
-4c, 

(52) 


(S') 

We are left with terms related to There is the finite contribution form the pion- 


nucleon vertex renormalization as discussed in section O!. It leads to as given 

in eq. (|3^ . The remaining contributions will be estimated by resonance exchange. This 


works well in the meson sector |^, it is, however, more complicated in the baryon case 
since one can have excited nucleon intermediate states as well as scalar and vector meson 
couplings to two or three pions. For our case, we estimate the genuine counter terms from 
single and double X* excitations and meson exchange, see Fig. 6 . Notice that 
at threshold, the vector meson contribution analogous to the scalar meson one vanishes. 
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To be precise, the p has a P-wave coupling and the Stto; vertex vanishes at threshold. 
In principle, there could also be a contribution due to the Oi Svr coupling. However, the 
branching ratio Oi —>■ (tttt)^tt is so small that we can safely neglect this contribution. We 
now turn to the non-vanishing parts. First, we consider the A(1232). If one treats it non- 
relativistically, all graphs with double-A excitations vanish at threshold, as it is e.g. the 
case in the model of ref.[^. These diagrams are obviously of order and proportional 
to (fj, which vanishes at threshold. Relativistically, the double-A graphs hrst contribute 
at order to D 12 . 

Second, there is the Roper, A*(1440). It has a very large width for decaying into a 
nucleon and two pions. The corresponding Lagrangian CN*NTnT for S-wave emission is 
discussed in detail in appendix ^ since in the available literature it is not treated in its 
most generality (as needed here). On the other hand, Cn^n-k is standard, 

jCn^n^ = 7 m 75 + h.c. , (53) 

with y/R = 0.53 ± 0.04 from the total width (calculated relativistically) .Q Putting 
pieces together, the Roper contribution is 


Df = 0 , = 


N* 


{cI + c^)9aVRM^ 






m] 


(54) 


The coupling constant combination c^ + C 2 = —1.56 ±3.35 GeV“^ follows from the partial 
width T{N* — N{'K'k)s) with both pions in an S-state [^. The details can be found in 
appendix p. 

The next resonance which decays into a nucleon and two pions is the Zli3(1520). Non- 
relativistically, this D-wave state cannot be excited by the initial P-wave ttN system and 
thus we neglect such a possible contribution. Of course, at energies above threshold such 
resonances become important, as witnessed e.g. by the study of the reaction yp —>• 7 r+ 7 r“p 
in [^. We remark that our treatment of the baryon resonance contributions is in good 
agreement with the partial wave-analysis of ttA —>■ tttiN by Manley et ah In their 
table VI one sees that at the lowest energy considered, only the Roper leads to a sizeable 
cross section. There is no sign of the A and the ZI 13 contribution is approximately one 
order of magnitude smaller than the one from the Pn. This lends further credit to our 
argument of exclusively keeping the Roper excitation to estimate the baryon resonance 
contribution to the threshold ttttA amplitudes. Third, there is scalar meson exchange. 
Its main contribution comes from diagrams with a nucleon pole. This is supposedly 
contained in the empirical values of the q, as discussed in some detail in Refs.. 
The remaining terms are of the form depicted in Fig. 6 c and include the vertex irSNN. 
In the absence of any empirical indication about the strength of this coupling constant 
and the lack of theoretical models thereof, we will set such terms to zero. It is conceivable 
that this is a good approximation since the bulk of scalar meson exchange is expected to 
be encoded in some of the terms proportional to the constants q. 

use here the width as obtained from the speed plot, not the model-dependent Breit-Wigner fits, 
Btot = 160 ± 40 MeV Q and as branching ratios BR(A^* —> Ntt) = 0.68 and BR(iV* —> N{TnT)s) = 


0.075 ±0.025 
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5.5 The master formula 


We have now assembled all pieces to give the chiral expansion of the threshold amplitudes 
Di^ 2 -i or equivalently, Aw and . 4 , 32 . The chiral corrections to the amplitudes Di ^2 follow 

||5|,|6|,P,|8|, 


by combining eqs. 


We have 


A = + 0{Ml) , * = 1,2, 


with 




Q'kN 

8 mF^ 




nyT = 


Q'kN 


8mF3 


3 +yA 


(55) 


( 66 ) 


and we do not write down again the various pieces contributing to the order corrections 
explicitely. Symbolically, they take the form 


= £)4:Born) ^ ^(2,loop) ^ ^ jj{2,Ci) ^ ^{2,GTR) ^ ^{2,TV*) 


* = 1 , 2 , 


(57) 


with ^ = 0, and we neglect the imaginary parts as discussed before. 

The LETs were hrst derived in ref.[^. With the help of eqs. (|^ , |^) as well as eq.(§), we 
can make explicit the tttt scattering lengths (i.e. give the constants d^j, cf. eq. (pd]) ). This 
determines the relation between the threshold ttuN amplitudes and the tttt scattering 
lengths one order beyond the improved representation, eq.(|l0l), and reads. 


Re Ao = dvr —^ + Ao + 0 {Ml) 


m IQnFF^ 

0/1 O /To dirN I '^v/fOS'/rTvAfTT 


+ A 


32F5 

VWqaM^ 


128F5 


+ 0(M3 


(58) 

(59) 


with 


An = 16 


(sr,-ir, + 3F)-^(?^ + i2al) + 


UFl 


m 


2Ci — 2 c 3 + C 4 + Co ) + 


TT^ 


{28gj 


16) In —-—h —- :r 9 A + 0 -‘^ 9 a + t) V^ln(2 + a/S) — 321 

A b 3 4 


/ F'^ \ 

+ 96( —(c^ + c^^fR —-h 5ii ) 

V m — m ) 


mA — m 
64A 


(60) 


/ \ 2F^ 64A / \ 

A 2 = 32(4i[-2iA -^ (1 + 12 gl) + -^ (-4ci + 4 c 3 + C 4 + Sc'a + 4d!, ) + 

V / vn^ m \ / 


TT^ 


{4g\ — 4) In h g~ H) + (^S'a ~ + V^) + 40 / 


(61) 


where the factors in front of Ao ,2 have been chosen such that the Ao ,2 are numerically of 
the order 0(1) and I is given in eq. (p7D . For the discussion of the respective numerical 
values, we abbreviate the O(M^) contributions as follows. 


An 2 — Anb + An°™ -r 2 w i-io,2 




A 0,2 — ^ 0,2 T =^ 0,2 

where the various terms can be read off form eqs.(^,R). Obviously, A^*^ = 0. 


(62) 
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6 Results and discussion 


First, we must fix parameters. Throughout, we use = 93 MeV, = 139.57 MeV, 
m = 938.27 MeV, = 13.4, gA = 1.26 and m* = 1440 MeV. This leads to bu = —0.012. 

Consider the low-energy constants q. From the formulae given in section |5.4| and 
with the empirical information contained in table 2.4.7.1 of together with (J^Ar(O) = 
45±9 MeV and the scattering volumina Cq = (0.208T0.003) M~^ and di = (—0.069± 
0.002) M~^ 1^, we arrive at the numbers given in table 1. We note that these numbers 
are typically a factor 1.5 smaller than the ones form the determination including the loop 
effects at order [M. The uncertainty in the table reflects the spread of the various 


determinations (to order if possible (like for ci, C 3 and C 4 ), otherwise the uncertainty 
of the empirical input. 


Cl 

4 

c" 

C 2 

C 3 

C 4 

-0.64 ±0.14 

-5.63 ±0.10 

7.41 ±0.10 

-3.90 ±0.09 

2.25 ±0.09 


Table 1: The low-energy constants q. All values in GeV ^ 


From these numbers, we deduce 


= 9.54 ± 0.21 , = 1.16 ±0.50 . 


(63) 


Next, we need the low-energy constants from the meson sector, f '1234 
from the recent analysis of 77^4 data beyond one loop and the £ 3^4 
paper [H, 


. We take the £ 1,2 
from the classical 


10^ 7^(lGeV) = -5.95 ± 1.06 , 10=^ 7^(lGeV) = 4.35 ± 2.75 , 


10^ £^(lGeV) = 1.64 ±3.80, 10^ £^(lGeV) = 2.29 ± 5.70 , (64) 

leading to 

Ao' = -0.93 ± 0.40 , A 2 ' = -1.04 ± 0.22 . (65) 

at the scale A = 1 GeV. If one lets A run from 0.5 to 1.5 GeV, one hnds Aq = —0.49... — 
1.05 and A 2 = —1.06 ... — 0.98. To acount for this, we scale up the uncertainty of Ag by 
a factor 1.5. The Born contribution to Ao ,2 is small, we hnd 


^Born _ _g_34 ^ ^Born _ _g_3g _ 

The corresponding loop contributions are also readily evaluated, 

AIOOP ^ ^ ^ ^loop ^ ^_g^ ^ g_^g _ 


for A = 1 GeV as the central value and the uncertainty accounts for the variation if A 
varies from 0.5 to 1.5 GeV. To end this part, we give the GTR and Roper contributions, 

= A^™ ± A^* = 0.22 ± 3.58 , ( 68 ) 
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where the large uncertainty stems from the Roper couplings c* + c^. From these numbers 
we can already conclude that the 1/m corrections to the ttN scattering amplitude (~ 
^ 0 * 2 ) ^ dominant role in the total correction of order together with the Roper 

excitation. Also, if one separates out the tttt interactions, their contribution is comparably 
small. This indicates that a very accurate determination of the S-wave ttyt scattering 
lengths will be very difficult (see below). 


Consider now the chiral expansion of the Hi, 2 - The LET values based on eqs. (|56|) are 
£^let _ 2A1 fm^ and = —6.76 fm^ compared to the empirical ones 
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HfP = 2.26 ± 0.12 fm^ = -9.05 ± 0.36 fm^ 


(69) 


These numbers result from a best £t to the near threshold cross section data for 


TT^Ti^n and tt' 


■p 


TT^Tl^n 


|Tl| IT^ using the correct flux and three-body phase space 


factors. The LET predictions show the expected pattern of deviation from the empirical 
values (compare the discussion in Refs. [pA|] |]I^). The various 0{M^) corrections to Hi ,2 
are summarized in table 2. 


dP 

Born 

Loop 

ii 

Ci 

GTR 

N* 

1 

-0.08 

0.25 ±0.13 

-0.17 ±0.18 

0.24 ±0.10 

0 

0 

2 

0.15 

0.09 ±0.48 

0.39 ±0.21 

-2.85 ±0.06 

0.32 

-0.40 ±0.90 


Table 2: Various corrections to Hi 2 - All values in fm^. 

TT 

Adding the uncertainties in quadrature, the predictions for Hi ,2 to order are 

Hf^ = 2.65 ± 0.24 fm=^ = -9.06 ± 1.05 fm^ . (70) 


These results for Hi ,2 are compatible with the empirical values, eq. (|6^ , within one stan¬ 
dard deviation. We notice that there are large cancellations in the contributions to 
Hi, whereas the chiral corrections to H 2 at this order are clearly dominated by the terms 
proportional to the low-energy constants Cj and the Roper excitation. In particular, one 
reads from table 2 that the loop contribution to the / = 0 amplitude ^10 is rather small, 
roughly -4% of the LET value. In contrast to this, the loop corrections for the / = 0 ttyt 
scattering length Oq are sizeable, about 25% of the current algebra value (at A = 1 GeV). 
This signals that contrary to expectations, the reaction ttN —>■ ttuN at threshold is not 
very sensitive to the tttt final state interactions. Also, if we disentangle the terms of order 
MH) 

{n = 0,1, 2) we hnd that the convergence of the expansion for Hi is good whereas 
in H 2 one still hnds sizeable corrections at u = 2, 


£>i = 1.59 • (1. + 0.52 + 0.15) fm’ , D 2 = -4.76 • (1. + 0.42 + 0.48) fm’ 


(71) 


for the mean values of Hi, 2 . This shows that the chiral expansion for Hi is converging. 
Matters are different for H 2 . One has at least to calculate the terms of order before 
one can draw a clear conclusion about the accuracy with which H 2 can be calculated. 
This, however, goes beyond the scope of the present paper. 


17 

























We now turn to the determination of the S-wave tttt scattering lengths. ^ Clearly, due 
to the large corrections and uncertainties (due to the Roper) in D 2 , deducing from 
the master formula can only give an indicative result. We hnd 


= 0.21 ±0.07 


(72) 


to be compared with the current algebra value of 0.16 and the CHPT prediction of 0.20 ± 
0.01 0. We notice that the theoretical prediction has a much smaller uncertainty than 
the number extracted from the titiN threshold amplitude. 

Matters are different for Oq since this quantity is entirely sensitive to Di. Adding the 
empirical and theoretical uncertainties in quadrature, we hnd 

a2 =-0.031 ±0.007 , (73) 


consistent (within one standard deviation) with the one-loop CHPT prediction of Gasser 


and Leutwyler, Oq = —0.042 ±0.008 [|| 0. We remark that for the combination 2ao 


■ Soq 


we have 0.577 ± 0.144, consistent with the universal curve, (2aQ — 5aQ)univ.curve = 0.614 ± 
0.028. This clearly indicates that previously determined scattering lengths based on the 
Olsson-Turner model [|T^ lITOIl should not be trusted. The lesson to be learned here is that 


even with an improved calculation there will remain sizeable theoretical uncertainties 
which will make a more accurate determination of the isospin zero S-wave tttt scattering 
length from the tttiN threshold amplitudes very difficult. In contrast, one can hope to 
sharpen the determination of Oq, eq 


7 Summary 

In this paper, we have considered the reaction ttN —>■ tittN at threshold in the framework 
of heavy baryon chiral perturbation theory. The pertinent results of this investigation can 
be summarized as follows: 


We have calculated the chiral expansion of the threshold amplitudes Di and D 2 (or, 
equivalently, to Aw and ^, 32 ) up-to-and-including the quadratic order in the pion 
mass. This amounts to the hrst corrections to the low-energy theorems derived 
in Ref.||17||. The resulting values for D 12 agree with the empirical ones within one 
standard deviation. For Di, the 0{M^) corrections are small, the corresponding 
corrections to D 2 are sizeable. The latter are mostly related to chiral corrections to 
the 7 tN amplitude and the excitation of the A^*(1440) resonance. 


• Based on this improved representation for the threshold ttttN amplitudes, one can 
deduce the isospin two S-wave tttt scattering length, Og = —0.031 ± 0.007. This 
number is compatible with the one-loop chiral perturbation theory results. However, 
the ensuing uncertainty is still sizeable and it will be difficult to further improve 

add the uncertainty from the theoretical determination and the one from the fit to the data in 
quadrature. 
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upon it. Due to the large corrections in D 2 , one can only deduce a broad range 
of values for Oq from this calculation, af] = 0.21 ± 0.07. The point here is that the 
threshold tttiN amplitudes are much less sensitive to the four-pion vertex than to 
other effects like resonance excitations and uncertainties in the ttN amplitudes. At 
present, it appears that the threshold n-nN data are best suited to pin down the 
isospin two, S-wave ttyt scattering length. 


As a by-product, we have found a new coupling for the A^*(1440) decay into the 
nucleon and two pions in the S-wave. It differs from the conventionally used one 
35|| through its explicit energy-dependence (i.e. factors of the pion energies), see 
appendix p. It would be important to disentangle these two couplings to reduce 
the uncertainty related to the Roper excitation. This could eventually be done 
by photo-exciting the Roper and study the decay N* N{ 711^)3 in the threshold 
region (i.e. for s/s > m*). 


Finally, we stress that an order calculation should be performed to further tighten 
the chiral predictions and to get a better handle on the chiral expansion of D 2 . Also, 
a consistent calculation to treat all isospin violating effects (like e.g. — m„, virtual 
photons and alike) has to be performed. 
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A Loop functions 


Here, we present explicit formulae for those loop functions entering the calculation of 
ttN — ttitN at threshold which are not given in the review [Q. All propagators are 
understood to have an inhnitesimal negative imaginary part in the denominator {—ie). 
We use dimensional regularization to compute divergent loop integrals and expand them 
around d = A space-time dimensions. In the reaction tiN —>■ irnN at threshold we 
encounter the following loop integrals involving three and four propagators, 




1 ) 


dH 


i J {27rY P) (M2 - (/ + 2g)2) 

1) 


— ^> 0 ( 0 ;), $ 3 ( 0 ;) -t-... (A.l) 


= A/o{u) , + 


(A. 2 ) 


i J {27 tY V ■ Iv ■ {I + 2q) (M^ — P) (M2 — (/ -|- 2^)2) 

with uj = V ■ q , and the ellipsis stands for terms proportional to q^qi,, ■ ■ ■ 

which are not needed in the actual calculation. For our purpose the loop functions $ 0 , 3 ( 1 ^) 
dehned above are to be evaluated at a; = ±M.,r + *0, 


MM.) = 




7T 


a/s ln (2 -|- a/s) -|- i ttVs 


M-M.) = 


L2 


TT 


+ y31n(2 + y3) 


(A.S) 

(A.4) 


^d,{M^) — 2 Mt^ L -|- 


itL[21„^+--A V31n(2 + 73) - t7rV3 
167r2L A 6 S ^ ^ 


$ 3 (—M^) — —2M„L + 


M. 

167r2 


-2 In ^ + ^ + I “ + Vs) 

A 0 S 


(A.S) 


(A. 6 ) 


with L dehned in eq.(|^. In order to get the values of the functions \I'o, 3 (<^) at cu = 
one simply uses the identity 


v-{l + 2q)v-l 2v-q\v-l v ■ {I + 2q) 
shifts the loop momentum and hnds 


MM.) = ^[MM.) - M-M.)] = - 21 n (2 + Vs)] 


(A.7) 


(A.S) 


^3(M^) = ^[<F3(M,)-<F3(-M^)] = 2L+ ^ 


2M^ 


167r2 . 


2 In —— - + a/s ln (2 -|- Vs) —i— a/S 

AS 2 

(A.9) 
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B The decay Ar*(1440) —^ N{7T7t)s 

In this appendix, we discuss in detail the Roper decay into the nucleon and two pions 
in the S-wave since that is not treated in generality in the present literature. In fact, to 
order the pertinent Lagrangian Cn^Ntt-k contains (at least) two terms. 


C-N*N-k-w — X+ — 3^2 ’ 


m 


(B.l) 


with '^n,n* the relativistic spin-1/2 helds. In fact, the second term is not unique, but for 
the non-relativistic formulation, all relativistically inequivalent forms lead to the same 
operator proportional to (n ■ m)^. Using 


X+ M^(2 p2+•■•)’ 


Uf, = + 

Utt 


(B.2) 


one hnds that the hrst (commonly used) coupling is energy-independent whereas the 
second depends on a;ia; 2 , the product of the energies of the two pions. Denoting by Fg 


the partial width r(A^* —> N{ 7111 ) 3 ), we hnd using eq. ([B.l|) , 
3 


Fg = 


IbTT^Fi 


[ [ duJiduj 2 {Tn* + m - UJI - UJ 2 ){c*M^ + cl u!iU! 2 )‘^ , (B.3) 

J Jz^<l 


with 


a;ia;2 — m*(a;i + U 2 ) + + (m*^ — m^)/2 


z = 


(a;f-M2)(n;i-M2) 
Performing the integrals in eq.( p.3|) , we have 

Fg = 


0.498cf + 2.708ctc2 + 3.714cf 


10"^ GeV^ , 


(B.4) 


(B.5) 


which is somewhat surprising since one would not expect the factor of pion energies 
leading to such a difference compared to the energy-independent coupling. Eq. (|B.5|) 
dehnes a rather elongated ellipse. The determination from the width does not £x the 


signs, we choose it to be negative for since the hrst term in the Lagrangian eq. (|B.l|) 
can be considered as ’’transition cx-term” (analogous to the Ci term in the dimension 
two pion-nucleon Lagrangian) and we know that Ci < 0. Similarly, we set C 2 > 0. For 
the threshold amplitude D 2 , only the sum of these two coupling constants is relevant. 
In the absence of any further empirical information, we take a mean value between the 
two extrema, = 0 or C 2 = 0, and allow for a large uncertainty to accomodate both 
possibilities. This leads to 


c* + c* = -1.56 ± 3.35 GeV^ , (B.6) 

which we use in the main text. These two couplings could be disentangled by a precise 
study of the energy dependence of the Roper decay in the threshold region. 
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C Calculation of -D 1 J 2 the heavy mass formulation 

In this appendix, we briefly sketch the derivation of the contribntion -0^25 6qs.(^|S), 
within the framework of the heavy nncleon Lagrangian. The starting point is the effective 
heavy nncleon action to order (for details, see appendix A of Ref.[pO[|), 


= d^xN 


2 m 


N , 


(C.l) 


where the ellipsis stands for terms not contribnting to the Ci/m corrections. The and 
are qnantities of order g*. While the A^^'^ connect the “large” components of the heavy 
nncleon fields, the give the transitions form the “large” to the “small” components 
before one integrates ont the latter (for details, see ref.|^). There are three different 
strnctnres contribnting to (from now on, we drop the snperscript ’cj’)- First, there 
are terms of the type 

(C.2) 
(C.3) 
(C.4) 


which lead to 

with 

Second, the combination 
contribntes to Di and D 2 , 


^(2) 1 ^ 

V ■ i 2 m 

Di = 2 P (—Cl + C2 + C2 + C3) , 
9aMI 


P = 


mP^ 


2 m 


^o^(i)Vi3(2)+h.c. , 


Di = -P (-C2 + C4) B2 = -P C4 


and third, there is a D 2 contribntion from 


^(3) 1 ^(1) ^ 
V ■ I 


which reads 


D2 = -P{c’2 + 24] 

Pntting pieces together, we arrive at eqs.(ff^,|l). 


(C.5) 

(C.6) 

(C.7) 

(C.8) 
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Figure Captions 


Fig.l Bom graphs. Solid and dashed lines denote nncleons and pions, respectively. The 
upper two diagrams are the contact and pion-pole graphs. The others are suppressed 
by powers of 1/m in HBCHPT. 

Fig.2 Renormalization of iG and to one loop. Crosses denote counter term inser¬ 
tions. The double and wiggly lines represent the pseudoscalar density and the axial 
current, in order. For other notations, see Fig. 1. 

Fig.3 Renormalization of the pion-nucleon vertex. Counter term insertions are not shown. 
For notations, see Fig. 1. 

Fig.4 One-loop diagrams non-vanishing at threshold. Graphs with the two out-going 
pion lines interchanged are not shown. For notations, see Fig. 1. 

Fig.5 1/m suppressed contributions from the relativistic dimension two pion-nucleon La- 

( 2 ) 

grangian. The circle-cross denotes an insertion from proprtional to one of the 
low-energy constants q. For notations, see Fig. 1. 

Fig.6 Resonance saturation, (a) and (b) refer to nucleon excitations, like the A(1232) 
or the iV*(1440) (double lines) and (c) to a t-channel meson exchange (double line 
denoted ’R’). 
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Figure 5 
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Figure 6 
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